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Motivated by the lack of an obvious spectroscopic probe to investigate non-conventional order such 
as quadrupolar orders in spin S > i systems, we present a theoretical approach to inelastic light 
scattering for spin-1 quantum magnets in the context of a two-band Hubbard model. In contrast 
to the S — \ case, where the only type of local excited state is a doubly occupied state of energy 
U, several local excited states with occupation up to 4 electrons are present. As a consequence, 
we show that two distinct resonating scattering regimes can be accessed depending on the incident 
photon energy. For hwi„ < U, the standard Loudon-Fleury operator remains the leading term of the 
expansion as in the spin-| case. For huui„ < 4U, a second resonant regime is found with a leading 
term that takes the form of a biquadratic coupling ~ (Si • Sj) 2 . Consequences for the Raman spectra 
of S=l magnets with magnetic or quadrupolar order are discussed. Raman scattering appears to be 
a powerful probe of quadrupolar order. 



I. INTRODUCTION 

The theoretical quest for exotic phases of magnetic 
matter over the last few decades has opened the way for 
a systematic investigation of models by various analyti- 
cal and numerical means and has led to the conceptual 
understanding of different states. This is for instance the 
case of resonating valence bond states which have been 
investigated in the context of high-Tc superconducto rs.^ 
frustrated spin systems^ and quantum dimer models! 3 -^ 
Yet, suggesting a technique to unravel the experimental 
fingerprints of an exotic phase often remains a challenge: 
for instance, in the spin liquid case, it is only recently 
that theoretical proposals have been made to detect this 
state on the kagome lattice with Raman scatteringSEl and 
that experiments have been carried outP 

Raman scattering was discovered in the beginning of 
the twentieth century. Loudon and Fleury showed how 
to use it to detect magnetic excitations in the beginning 
of the sixtiesPlt was used for instance in the late eighties 
to estimate the exchange in the cupratesPmore recently 
the investigation of magnetic properties with this spec- 
troscop ic techn ique in various contexts like frustrated 
systems ^ 6 ! 10 ! 11 ! or iron pnictideJ^ has been a very ac- 
tive field of research. In fact, inelastic light scattering 
takes advantage of the fact that the response is in essence 
linked to many-body physics and that the photon polar- 
ization enables t o coll ect valuable insights on correlated 
electron systems! 13 * 14 ^ Following this route, we present 
here inelastic light scattering as a natural probe to char- 
acterize the order and the low-energy excitations in spin 
5=1 models. 

The pecularity of quantum magnets with 5 > 1/2 is 
that they can break SU(2) symmetry with a local order 
parameter which is non magnetic. The simplest exam- 
ple is provided by quadrupolar order for 5 = 1 systems, 
where the local order parameter is not the spin but a 
rank-2 tensor. This has been extensively studied in the 



context of the bilinear-biquadratic model defined by the 
Hamiltonian: 
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Ferro- and antiferroquadrupolar phases have been i den- 
tified both on the triangula r! 15 * 16 ! and square^E^ lat- 
tices for sufficiently large biquadratic interactions, and 
the compound NiGa 2 S4, where Ni 2+ ions form a trian- 
gular lattice of spins 1, has been sugge sted to exhibit 
some kind of quadrupolar ordeA 5 -^ 6 * 19 * 20 !. Possible mech- 
anisms to produce large biquadratic interactions include 
spin-lattice coupling^ and situations with quasi-orbital 
degeneracy^. 

A direct observation of quadrupolar order remains 
a challenge however.^ Quadrupolar states being non- 
magnetic, conventional experimental techniques such as 
neutron scattering are insensitive to quadrupolar order. 
By contrast, we show in the present work that, since light 
naturally couples to the charge, inelastic light scattering 
offers an alternative to investigate non-magnetic states. 

The paper is organized as follows: in a first section we 
present the microscopic model and the derivation of the 
effective magnetic light-scattering operator. We discuss 
the form of this effective 5=1 operator compared to the 
more conventional 5 = \ case. We show that two dif- 
ferent resonant regimes are accessible, depending on the 
incoming photon energy and we discuss the different po- 
larization geometries relevant to our problematic. This 
leads to a section in which the expected spectra for the 
different phases are displayed, followed by a section de- 
voted to a discussion concerning potential experimental 
checks, and by a brief conclusion. 
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II. EFFECTIVE LIGHT SCATTERING 
OPERATOR 

A. Microscopic parameters 

From the point of view of purely atomic physics, spin- 
1 states can be achieved in the case of transition metal 
ions, for instance Ni 2+ (3d 8 ) in a cubic environment that 
leads to the standard t2 g — e g splitting of the orbital d- 
shell. The intra-atomic electron-electron interaction of 
the partially filled shell can then be reformulated as a 
Hund's coupling which favors states maximizing the spin. 

As we want to keep a very general Hubbard model, we 
will consider here a lattice of such sites: two degener- 
ate orbitals at half-filling, i.e. 2 electrons per site, with 
nearest-neighbor hopping and on-site interactions that 
include inter- and intra-orbital coupling. This two-band 
Hubbard model is described by the following Hamilto- 
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where i,j are the site indices, m, m! refer to the orbitals 
and a to the electronic spin. The hopping parameters 
between two neighboring orbitals is t m m r, the on-site 
Coulomb repulsion is denoted by U rnm ' , Jh represents 
the Hund's coupling and Jp the pair hopping amplitude. 
For the sake of clarity, we restrict ourselves here to a 
square lattice, although the extension to other lattices is 
straightforward. Furthermore, we assume that the addi- 
tional relations, typical of cubic symmetry, are satisfied, 
namely: U aa = U bb , U = U aa - 2J H and J H = 2J P . 

In the Mott insulator regime, at half-filling and for 
U m m' 3> tm.m' , the electrons are localized and Hund's 
coupling favors triplet states on each site. The relevant 
degree of freedom is then a spin S = 1. At second or- 
der in perturbation theory, the resulting effective spin-1 
Hamiltonian is the standard Heisenberg model: 
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As in the case of S = 1/2 effective models for single-band 
HubbardJ 23 * 24 ! pushing the perturbation to fourth order 
leads to the emergence of additional terms such as 4-site 
terms 25 and biquadratic interactions. 



B. Derivation of the effective operator 

1. General method 

Inelastic light scattering techniques like Raman are 
photon-in photon-out techniques where the incident pho- 
ton couples to the charge. Shastry and Shraiman gave a 
microscopic description of the process in Ref. [26 The 
procedure is however given in a more pedagogical way 
in Ref. [B] As we follow very closely the latter, we will 
not explain in great detail the derivation but just sketch 
the main ideas. The photon-electron interaction is intro- 
duced via the Peierls substitution: 



-> Ci, 



■ exp 
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where A is the photon vector potential. 

For incoming photon wave-lengths much larger than 
the lattice spacing, it can easily be shown that the cou- 
pling of the microscopic Hamiltonian to the photon af- 
fects only the kinetic terms and generates a current that 
depends on the incident (ei n ) and scattered (e ou t) po- 
larizations. After second quantization of A, a and a' 
being respectively the photon creation and annihilation 
operators, the current term reads: 
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V, where V stands for the 
j is the vector connecting site i and site 



7 , where H-y represents the 



with g m ,out = \/hc 2 /uj^ 
volume, and e 
j- 

Hence, collecting all the parts of the Hamiltonian de- 
scribing the electronic system coupled to the light we 
have: % = Hub + He + H 
purely photonic part of the Hamiltonian 

In the limit U ^> t, the electrons are localized and the 
low-energy spectrum of the system can be described by a 
spin Hamiltonian. Furthermore, in the limit \U — LOi n \ < 
t, H c can be treated as a perturbation and one can de- 
rive an effective magnetic Raman operator in the spin 
sector. Yet, if one is interested to have an effective oper- 
ator beyond the standard limit \U — Wi n \ <§C t, one has to 
push the perturbation theory to fourth order. It is worth 
pointing out that this effective description conserves the 
quantum numbers of the original Hamiltonian and scat- 
tering operator: for instance since Hub and He conserve 
the total spin, so will the effective Raman operator and 
thus we recover the fact that AS = for Raman excita- 
tions. 

The effective scattering Hamiltonian obviously de- 
pends on the chosen polarizations of the incoming and 
scattered photons; to clarify the discussion we present in 
FigjTJthe different geometries we have investigated. 
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Figure 1: Polarization geometries for the incident and scat- 
tered photons. The present decomposition refers to the scat- 
tering sectors at second order in the derivation. Higher order 
terms will mix the sectors as indicated in Table I of Ref. 1261 



The scattering operator can in general be decomposed 
in two different ways: the first one involves the symme- 
try of this operator and the second refers to the differ- 
ent incoming and outgoing polarization vectors. These 
two decompositions are equivalent at second order, but 
they are different at higher order (see Ref. I2"?)l . Here, we 
choose to decompose the scattering operator with respect 
to the polarization geometries shown in Figjl] To stick 
to the most conventional notations, the top-left crossed 
polarization geometry of Fig (T] will be referred to as Bi g , 
while the bottom- left geometry will be called A\ g . 



2. Single-band case: results for 5 = 1/2 



This calculation has already been done in Ref. [2"6"1 
However, Ko et aZP recently pointed out some differ- 
ences in some channels between the operator calculated 
initially and their own calculation. So, as a warm up, 
we have rederived the operators for every channel. We 
found further differen ces at fourth order with respect to 
the initial calculation. 1221221 Denning V a = J2 r S r ' S r+a , 

Qa,P,S = J2r( S r- S r+ a )(S r +l3-S r+S ) and A = t/(U-U)i), 

we found: 



Blg = 4tA ( - + 2A ) (V y - V x ) ~ 2tA A {V 2y - V 2xj 



Alg = UA ^- + 6A 2 j (V y + V x ) - 2tA A (T 2y + V 2x ) 
- 8tA 3 {V x+y +V x -y) 

"I" 32tZ\ (.Qx,y,x-\-y ~t~ Qy.x.,x-\-y Qx-{-y,x,y^ 
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A ' lg = AtA(l + 8A 2 )V x + 16tA 3 V v 
- AtA 3 V 2x - 8tA 3 {V x+y + V x -y) 

+ 32tA (Qx,y,x-\-y H~ Qy..x,x-\-y ~~ Qx-\-y,x,y) 
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In Appendix|X] we discuss the differences between this 
operator and the result of Ref. 26 28 after rewriting it 
following the symmetry-based decomposition used in this 
reference. To test the validity of our operators, we have 
also compared on small clusters the Raman spectra of the 
original single-band Hubbard model with those obtained 
with the effective operators, with the conclusion that the 
spectra obtained with our operators agree much better. 



3. Two-band case: results for S — 1 

The general procedure is similar to the one of the 
single-band case, but starting from the two-band Hamil- 
tonian described in Subsection III At There are two main 
differences between this case and the single-band case. 

First, in the one-band case, the only excited states that 
were considered were the terms with two particles on one 
site, inducing a Coulomb repulsion U. These terms will 
still be present in the two-band case, but other terms 
where only the Hund's coupling is not satisfied will also 
be present. At fourth order, this leads to resonances at 
different incoming energies. 

Moreover, since spin-1 live in a 3-dimensional space, 
there arc eight non-trivial hermitian operators acting on- 
site (instead of three in the spin- 1/2 case). We thus ex- 
pect new types of operator to appear, as for example the 
biquadratic coupling {Si ■ Sj) 2 . 

To express the Raman operator in this case, we take 
the same definition for V a and Q a ,/3,8 as for spin-1/2 and 
we introduce in addtion the operators lZ a = y^ r (S r ■ 

Sr+a) 2 and T a ,/3 = J^ri^r ' S r + a )(S r + a ■ Sr+p). To 

fourth order in perturbation theory, the resulting opera- 
tors for Big and A\ g read: 
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+ A d (V X + y+V X -y) 

+ A h2 {V 2x + Vly) + A 3 (%,2x + T y .2y) 

4~ A^d (T x . x -\-y "F Tx } X—y "F 7y, x +y T—y^ x — y) 

~f" Ap (Qx ,x+y ,y ~F Qy,x,x+y Q^+3/,x,y) (4) 

The different coefficients are functions of the different 
parameters of the initial Hamiltonian, Eq.|2]). They are 
given in Appendix [B] 

C. Relevant limits and geometries 

One of the major advantages of Raman scattering lies 
in the possibility of using light polarization to select and 
characterize the excitation that one is willing to investi- 
gate. As the derivation remains very systematic, so far, 
we have taken care of all the different polarizations. At 
second order, only A\ g and B\ g geometries have non- 
vanishing operators. Therefore, in the rest of the paper 
we will mainly focus our attention on these two cases. 

In addition to the usual decomposition due to the in- 
coming and outgoing polarizations, one should also be 
aware of the information that can be accessed through a 
suitable choice of the incoming photon energy. The Ra- 
man operators and the associated prefactors of Appendix 
[B]lead to many terms. However, since all the prefactors 
depend on the incoming photon energy, it is possible to 
adjust u!i n to get close to a resonance and highlight spe- 
cific terms of the scattering operators. 

In the limit hu>i n < U + 2Jjj ~ U, the second order 
processes will be dominant, and the other terms can be 
neglected. Therefore, in this case, we can restrict our- 
selves to the Fleury-Loudon Raman operator: 



@Alg OC (S r • S r + X + S r ■ S r + y ) 

S r -\- X S r ' S r -\-y ) 



Another interesting limit occurs when huji n ~ AU . In 
this case, we will favor processes where there is an in- 
termediate state with four electrons at the same site. 
Such processes occur at fourth order if they involve only 
two sites. Moreover, as they can change the local spin 
by AS = 2, they have to be related to the operator 
(Si ■ Sj) 2 . This explanation is confirmed by the exact 
coefficient of the Raman operator given in Appendix [B] 
So in this limit, we can consider that the Raman operator 
reduces to: 
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Figure 2: Raman spectra of the spin-1 Heisenberg model 
(6 — of the Hamiltonian of Eq. [T]) for A\ g and B± g po- 
larizations. The left column represents the inelastic response 
as a function of the energy loss for an incoming photon of 
energy hu in < U, which corresponds to the Loudon-Fleury 
operators Ofl oc £\ [Si ■ S i+X ± Si ■ S i+y }. The right col- 
umn corresponds to an incoming photon energy of hujin < 4J7, 
and hence to the biquadratic form of the scattering operators 
Osiq oc J2i [( s i ' Si+x) 2 ± (Si ■ S i+y ) 2 ] . The spectra have 
been computed for a 16-site cluster. 



So, in two different limits, the Raman operator can 
be written in a quite simple way for the Bi g and the 
Aig geometries. The next section will be devoted to an 
investigation of the Raman spectra calculated from these 
four operators in different situations. 



III. RESULTS: EXPECTED SPECTRA 

The goal of this section is to look at the Raman re- 
sponse for specific angles 9 of Eq.Q. We do not mean 
here to give a quantitative and exhaustive analysis of the 
expected spectra but rather a qualitative description of 
them and, most importantly of their evolution through- 
out the different phases. These investigations are carried 
out numerically. The presented spectra have been ob- 
tained with Lanczos and the continued fractional for 16- 
site square lattice cluster with periodic boundary condi- 
tions. This cluster has additional symmetries^!, but this 
is of no consequence for our present purpose, which is to 
investigate the overall shape and behavior of the spectra 
for different phases. 

The results and spectra discussed in the present section 
are summarized in Fig(2]and in Fig|3] 



A. Heisenberg AF phase 

The most conventional Heisenberg antiferromagnctic 
phase is obviously reached for 8 — in the Hamiltonian 
of Eq.Q. The Raman response in the B\ g geometry for 
the Fleury-Loudon operator exhibits a bi-magnon peak 
which can be understood with a standard spin-wave cal- 
culation: Linear spin-wave theory leads to a peak at a 
Raman loss Ml = h (uji n — uj out ) ~ 8 J which softens 
down to hfl ~ 7.8J after taking into account magnon- 
magnon interactions, as shown for instance by Chubukov 
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Figure 3: Raman spectra for a model with purely biquadratic 
coupling (6 = — ^ of the Hamiltonian of Eq. Un for Ai g and 
Bi g polarizations. The left column represents tne inelastic re- 
sponse as a function of the energy loss for an incoming photon 
of energy fuo in < U , which corresponds to the Loudon-Fleury 
operators Ofl oc Y\ [Si ■ S i+X ± Si • S i+y }. The right col- 
umn corresponds to an incoming photon energy of TuuJin < 4[7, 
and hence to the biquadratic form of the scattering operators 
Osiq oc J2i [( s i ' Si+x) 2 ± (Si ■ S i+y ) 2 ] . The spectra have 
been computed for a 16-site cluster. 

and FrenkeP^ or Canali and GirvinP^l This is in qualita- 
tive agreement with our numerical result. For the same 
polarization, the biquadratic operator gives the same re- 
sponse with additional peaks for higher Raman losses 
corresponding to multi-magnon excitations. 

In the Ai g geometry, the Fleury-Loudon operator com- 
mutes with the Hciscnbcrg Hamiltonian and the Raman 
response vanishes. This is not the case for the biquadratic 
operator: as in the Fleury-Loudon case, the two-magnon 
features disappear, but the multi-magnon peaks, which 
arise from the non-commutating part of the operator with 
the Hamiltonian, remain visible. It would be really in- 
teresting to check this prediction experimentally, namely 
that intensity shows up at larger Raman loss energy when 
the incoming light frequency is in the range AU jh. 

B. Ferro-Quadrupolar order 

For quadrupolar order, the relevant order parameter, 
Q (see Ref . , and its descendant observables like Q • Q 
can be written in terms of spin operators: 

Q, Qj = 2 (Si • S,) 2 + S, • Sj + Cst. (5) 

An easy point to investigate corresponds to 9 = —tt/2 
where the Hamiltonian reduces to its biquadratic part. 

As already discussed, there are two resonant regimes: 
the first one for an incoming photon energy huji n < U 
leading to the bilinear Fleury-Loudon operator, and a 
second one for Mj in < AU at which the leading term of the 
scattering operator is proportional to (Si • Sj) . Since in 
Eq.Q both bilinear and biquadratic couplings appear 
in the quadrupole-quadrupole observable, both resonant 
regimes offer the possibility to investigate quadrupolar 
order. 

Concerning the physics occuring at fej jn < U, the 
Fleury-Loudon scattering Hamiltonian is expected to 



give access to quadrupolar excitations, and this turns out 
to be the case: The Raman response exhibits a marked 
peak for a Raman loss of Ml ss 5 J, which can be under- 
stood as coming from a pair of quadrupolar excitations. 
Indeed, a flavor- wave^212il ca j cu lation gives the dispersion 
of the quadrupolar excitations across the Brillouin zone, 
as displayed in Figj4] with a 2 J excitation at (ir, 0), lead- 
ing to a peak close to a Raman loss of A3 . This predic- 
tion corresponds to the excitation of two non-interacting 
quadrupolar excitations. However, this simple picture 
does not take into account particle-particle interactions, 
which in the present case appear to shift the peak po- 
sition at higher energy. At first glance, it is not obvi- 
ous that the Loudon-Fleury scattering operator offers the 
possibility to investigate quadrupolar excitations since 
this scattering resonance is mostly known to highlight 
the bi-magnon peak. Yet, in contrast to the conventional 
magnetic scattering, the two scattering geometries, A± g 
and Big , have comparable spectra, as none of the oper- 
ator commutes with the Hamiltonian. The similarity of 
the two spectra for the two polarization geometries seems 
to be a signature of a ferro-quadrupolar phase. 

The second resonance, for an incoming photon en- 
ergy of hujin < AU , is associated to a scattering op- 
erator with biquadratic spin terms. As quadrupolar or- 
dered phases are ground state of the Hamiltonian when 
(Si • Sj) 2 dominates, one expects that the effective Ra- 
man operator for hwi n < AU probes directly quadrupo- 
lar excitations. Indeed, this is confirmed by our calcu- 
lations and the obtained results can be explained in a 
straightforward manner. For the Ai g geometry, the Ra- 
man operator commutes with the Hamiltonian, hence the 
response vanishes. One should note however that, for real 
systems many terms appear in the scattering operator 
and longer-range terms in the Hamiltonian should also 
be considered; all this facts combine to give a finite Ra- 
man response even in this geometry. Yet, exactly like for 
conventional magnetic Raman scattering, the amplitude 
of the response should be much stronger for Bi g polariza- 
tions, as shown in Fig(2j One recovers the quadrupolar 
excitation at Ml 5 J. 



IV. DISCUSSION 

As summarized in Figs (2] & [3] the signal correspond- 
ing to a Neel ordered and a ferro-quadrupolar phase are 
qualitatively different. Hence, Raman scattering offers 
the opportunity to detect quadrupolar ordering in a rela- 
tively straightforward manner. The systematic is rather 
simple: one should analyze the Raman spectra obtained 
for two scattering geometries (parallel and crossed polar- 
izations) as presented here and tune the incoming photon 
energy, this last step being of course the most crucial one 
as it enables one to switch from the standard Loudon- 
Fleury limit to the biquadratic form of the scattering 
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Figure 4: Flavor-wave dispersion for a ferro-quadrupolar or- 
der corresponding to 8 — —ir/2 in Eq. |TJ. 

operator. This is one of the main messages of the paper: 
the microscopic derivation of an effective magnetic scat- 
tering operator for S — | systems leads to one resonance 
at huji n < U since there is only one possible interme- 
diate state. This situation is no longer valid for 5 = 1 
compounds, as a spin-1 is formed by a system of two elec- 
trons per site strongly coupled via Hund's rule, different 
intermediate states can be accessed: one with 3 electrons 
at one site, leading to the usual resonance at ftwi n < U; 
and another one with 4 electrons at one site leading to a 
second resonance at huii n < 4U. Of course, many other 
intermediate states are possible, however, at second or- 
der, the only possibility is to have 3 electrons at one site, 
therefore the Loudon-Fleury term will dominate; and, at 
fourth order biquadratic terms are only appearing along 
with intermediate states having quadruple occupancies. 

However, going from the first resonant regime (hu> in < 
U) to the second one (hu>i n < 4U) requires an adequacy 
between the materials and the available light source. In- 
deed, having a handle on the in and out polarizations re- 
mains much easier for visible light, and if U is too large, 
one would end up with incoming photons in the UV re- 
gion of the light spectrum. A promising route might be to 
try organic systems, in which all interactions, including 
U, are smaller than in oxides. 

Turning to the specific case of the ferro-quadrupolar 
phase, the experimental investigation would require two 
steps: i) for an incoming photon energy of Hu)i n < U, 
one should start with a geometry corresponding to the 
A\ g polarization and observe a spectrum with a peak 
at about Ml ss 5J. By keeping the same photon en- 
ergy and slowly rotating the outgoing polarization until 
getting to the Bi g polarization, one should not observe 
much variations but a small hardening of the main peak, 
ii) The second step consists in tuning the incoming pho- 
ton energy to hcj in < AU and to collect the spectra for 
the different polarizations from B\ g to A\ gi the Raman 
response should this time exhibit strong modifications: 
going from a marked peak to much broader features. In 
principle, performing the experiment within the Loudon- 



Fleury limit should be sufficient. However, if one does 
not know precisely the value of J in the considered com- 
pound, this two-step procedure allows a clear identifica- 
tion of the ferro-quadrupolar order. Precise calculations 
for the specific compound would also give information 
about J after fitting the data. Also, and from a more 
fundamental perspective, one should notice that the sec- 
ond resonance regime should exist even in the case of a 
standard Heisenberg phase. In this situation, the Raman 
spectra should exhibit not only a bi-magnon peak, but 
also more spectral weight at higher energy loss. 

V. CONCLUSION 

In the present paper, we have derived a general ef- 
fective inelastic light-scattering operator for spin-1. We 
have shown that this operator offers two different res- 
onant regimes depending on the choice of the incom- 
ing photon energy. On the basis of calculated spectra, 
obtained by exact-diagonalization of finite clusters, we 
have shown that the different phases of interest (Neel 
and ferro-quadrupolar) exhibit characteristic fingerprints 
that allow a clear identification of each type of ordering. 
This work, for the square lattice, illustrates the potential 
of Raman scattering technique as a probe for character- 
izing quadrupolar order; one of the best candidate for 
such a phase remains NiGa2S4 where the spin S = 1 lo- 
cated on the Ni 2+ ions form a triangular lattice. It is 
clear that a direct application for this compound is not 
entirely possible as it requires another derivation and a 
suitable decomposition of the scattering channels, which 
is left for further investigation. 
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Appendix A: Comparison of effective Raman 
operators 

In this Appendix we aim at comparing an effective 
magnetic Raman scattering operator in the Mott insu- 
lating state for a spin S = o model to the original Ra- 
man current operator that is associated to the Hubbard 
model. We present here a comparison for the Hubbard 
model on an 4 x 2 = 8-site ladder for B\ g polarization as 
well as for a 6-site chain. The effective scattering opera- 
tor derived in the present work showing some differences 
compared to Shastry and Shraiman's, we also provide 
graphs to compare our results to theirs. 

In order to stick to Shastry and Shraiman's original 
notations of Ref. 26, we remind their results to the reader: 
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O s = 2tA (V y + V x ) + 8iA 3 {V 2x + V 2y + V x+y + V x - y ) + 32tA 3 (Q x , y , x+y + Q y ,x 7 x-\-y Qx-\-y,x,y) 
O d = AtA [i - 4A 2 ] {V y - V x ) + 8iA 3 {V 2y - V 2x ) 

(Al) 

= 64iA 3 £ r ze^S r • (S r+A1 x S r+M /) 

O e = -mtA 3 {T x+y +T y - x ) 
While our derivation gives: 

O s = {2tA + 2AtA 3 ){V y +V x )-2tA 3 {V 2x +V 2y )-UA 3 {V x+y +V x - y ) 

-\- 32£A (,Qx,y,x-\-y ~l~ Qy ,x,x-\-y Qx-\-y,x : y) 

O d = AtA [1 + 2A 2 ] {V y - V x ) - 2tA 3 {V 2y - V 2x ) , 



Oo = 
O e = 

Ko et al) 2 * already noticed that there is no chiral 
term appearing in 01 unlike in Shastry and Shrainian's 
derivation. To check our derivation we perform a numer- 
ical analysis on finite clusters. The quantities we com- 
pare are the height of a most prominent peak in the exact 
model and the height of the same peak in the effective 
model at second and fourth order. We performed this 
analysis on two different clusters: 

1. A2x4 = 8-site cluster with open boundary condi- 
tions. This choice is motivated by the fact that it 
is the smallest cluster where some bounds have the 
correct prefactor in the x direction (for smaller clus- 
ters, some fourth-order renormalization factors will 
no be present in the term Si ■ Si+ X ) , while periodic 
boundary conditions would lead to the renormal- 
ization of some terms with respect to the infinite 
case coming from processes where a doublon goes 
through the four sites of one leg. The results for 
this cluster are presented in Fig. [5] 

2. A 6-site chain with periodic boundary conditions. 
This cluster is very usefull to check the prefactor of 
Si ■ Si+ 2x . As the term Si ■ Si+ X commutes with 
the Hamiltonian on this cluster, only the Si ■ Si +2x 
term contributes to the scattering amplitude, thus 
we can check very precisely the coefficient. Here, it 
is not possible for a doublon to travel through the 
entire system at fourth order, and therefore we can 
choose periodic boundary conditions. The results 
for this cluster are presented in Fig. [6] 

As can be seen in these figures, in both cases the spec- 



I 

trum derived from our effective Raman operator agrees 
very well with that calculated directly with the original 
Hubbard model, while the effective operator of Shastry 
and Shraiman leads to significant differences beyond the 
second order. 

1 1 1 1 1 1 

Effective model, second order 




0.99 0.991 0.992 0.993 0.994 0.995 0.996 0.997 



Figure 5: Comparison of the height of the main Raman peak 
for a fixed Raman loss as a function of the incoming photon 
energy for different models for a 2 x 4-site ladder with open 
boundary conditions. The hopping amplitude has been set to 
t = 5Hnn so that * < U and J = i£ = 



Appendix B: Coefficients of the effective spin-1 
Raman operator 



In this Appendix, we list the expression of the prefactors of the effective spin-1 Raman operator defined in Eq.(|4] 
The terms enclosed in a square box are those that dominate at the two resonances. 



Effective model, second order 
Effective model, fourth order, present work 
Effective model, fourth order, Shastry & Shraiman 
Hubbard model 




0.993 0.9935 0.994 0.9945 0.995 0.9955 0.996 0.9965 0.997 

co in /U 



Figure 6: Same as Fig. [5] for a 6-site chain with periodic 
boundary conditions. 
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where the parameters fi, £2, ^2p, ^4 are defined in term of the original microscopic hopping parameters by: 
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